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Defects of graphene are the most important concern for the successful applications of graphene since they 
affect device performance significantly. However, once the graphene is integrated in the device structures, 
the quality of graphene and surrounding environment could only be assessed using indirect information 
such as hysteresis, mobility and drive current. Here we develop a discharge current analysis method to 
measure the quality of graphene integrated in a field effect transistor structure by analyzing the discharge 
current and examine its validity using various device structures. The density of charging sites affecting the 
performance of graphene field effect transistor obtained using the discharge current analysis method was on 
the order of 10 14 /cm 2 , which closely correlates with the intensity ratio of the D to G bands in Raman 
spectroscopy. The graphene FETs fabricated on poly(ethylene naphthalate) (PEN) are found to have a lower 
density of charging sites than those on Si0 2 /Si substrate, mainly due to reduced interfacial interaction 
between the graphene and the PEN. This method can be an indispensable means to improve the stability of 
devices using a graphene as it provides an accurate and quantitative way to define the quality of graphene 
after the device fabrication. 



The electrical and physical properties of graphene have been extensively investigated for diverse applications 
such as electronic switch, sensors, transparent electrodes, fuel cells, and catalysts 17 . One of the key common 
challenges in these diverse applications of graphene is a wide variation in the material quality itself and the 
influences of the external environment via defects and surface reactions 8-12 . Thus, the origins of the instability of 
graphene devices and their sensitivity to external environmental factors have been extensively investigated 13-15 . 

For example, the conductivity of graphene is found to be greatly affected by many factors such as metallic 
residues originating from the growth processes, substrate materials, capping dielectrics, device fabrication pro- 
cesses such as patterning and cleaning, initial defect density, contact metals, operation ambient, temperature, and 
more. In addition, the operation history- dependent device drifts such as hysteresis, charge trapping, and inter- 
facial reactions also affect the characteristics of a device using a graphene. 

Since the defects in graphene and the influences of the environmental factors play a critical role in device-level 
operation 16 , methods to quantitatively monitor these factors, especially after device fabrication, are of utmost 
importance. However, macroscopic defects in graphene substrates have been primarily characterized using 
physical analysis methods such as visual inspections, Raman spectroscopy, atomic force microscopy (AFM), 
transmission electron microscopy (TEM), and scanning tunneling microscopy (STM) 17 ' 19 . The distance between 
defects, grain size, and relative density of sp 3 type defects can be measured using the shift in the relative position 
and width of the G-peak, D-peak, and 2D-peak of the Raman spectrum. Various TEM and STM analyses of grain 
size, defects, and contaminants have been reported in the literature 19-20 . 

While Raman spectroscopy, AFM, TEM, and STM are useful means in the analysis of the initial quality of 
graphene, these methods are not easy to use when examining graphene integrated in a device structure 21 ' 23 . For 
example, Raman spectroscopy becomes inaccurate when the graphene has a dielectric passivation layer. The 
beam size of a Raman spectrometer, which is typically larger than 0.5 urn, is too big for applications requiring 
small graphene patterns. Furthermore, Raman analysis cannot be performed with a thick metal electrode on 
graphene. AFM, TEM, and STM have similar limitations as they are destructive analysis methods. Thus, it has not 
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Figure 1 | A discharge current analysis method for graphene devices, (a), (b), Schematic showing the electrical connection and charge flow direction to 
measure the interfacial defect density for silicon MOSFETs (a) and the charging site density for graphene FET(b). (c), Schematic showing the 
sequence of discharge current generation in a single pulse cycle. Blue circle represents the intrinsic bulk carriers and red circle represents diverse defect 
sites, which has much slower charging and discharging time than the intrinsic bulk carriers. Step 1 is the initial state before applying a pulse. Holes are 
populated by the pulse (Step 2), Then, defect sites are charged up in step 3. When the pulse is turned off, majority carriers are discharged quickly(Step 4), 
and then trapped charge are slowly discharged(Step 5). Finally, channel charge density returns to the initial state, (d), I-V curve of a graphene FET 
showing the range of pulse bias. The discharge current is generated by a continuous pulse with a pulse height of AV. 



been possible to quantitatively analyze the quality of graphene at the 
device level, which is critically important for practical applications. 

Thus, the characteristics of graphene integrated in a device struc- 
ture have been primarily investigated indirectly using the electrical 
characteristics of devices such as current-voltage (I-V) or capacit- 
ance-voltage (C-V) characteristics 22 24 , which provide information 
on the mobility, location of the Dirac point, and hysteresis of the I- 
V curves 26 . Unfortunately, these device characteristics are known to 
vary significantly because they are very sensitive to environmental 
factors (substrate, capping dielectric, ambient, temperature, contact 
resistance, etc.) 27,28 . Since these aforementioned methods do not pro- 
vide a direct means of analyzing the quality of graphene after device 
fabrication, it is imperative to develop a new method to quantitatively 
identify the influence of defects and environmental factors even after 
device fabrication. Without a direct defect assessment method, it will 
be impossible to identify whether the changes in the device charac- 
teristics are induced by changes in the quality of the graphene or by 
the influence of extrinsic factors. 

In this work, a new electrical characterization method to obtain 
quantitative information on defects and other environmental factors 
representing the quality of graphene FETs has been proposed and its 
validity has been examined using devices fabricated on a large area 
monolayer graphene sheets with different levels of physical defect 



density (see Methods). This method, called the discharge current 
analysis (DCA) method, modulates the carrier concentration of gra- 
phene periodically using an external pulse bias and translates the 
frequency dependence of the charges discharged from the graphene 
channel into a density of charging sites. The validity of this DCA 
method is examined by correlating the electrically measured density 
of charging sites with the initial physical defect density of graphene 
measured with Raman spectroscopy. 

For silicon MOSFETs, a charge pumping (CP) method has been 
used to accurately analyze the interfacial defect density, which affects 
the device performance. The CP method fills up defect sites at the 
silicon-dielectric interface using minority carriers supplied from the 
source and drain at an inversion state, as shown in Fig. la. Then, 
majority carriers from the substrate are supplied to the interface by 
changing the silicon surface to an accumulation state. The majority 
carriers are supplied to the defects only in the channel because they 
cannot flow into the source/drain side due to the energy barrier of the 
pn junction. Finally, the defect density can be calculated by counting 
the number of majority carriers recombined with the minority car- 
riers in the defect sites. The measurement principle of CP is explained 
in detail in the Supplemental Section SI. Unfortunately, this method 
cannot be used for graphene FETs because there is no body contact to 
supply the majority carriers to the channel. As a detour, the carriers 
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charged and released from defect sites near the surface of graphene or 
graphene itself can be collected through the source/drain as shown in 
Fig. lb. This approach is feasible for a graphene because there is no 
barrier for the carriers due to Klein tunneling. However, the separa- 
tion of carriers from the defect sites and the channel itself becomes a 
technical challenge because majority carriers in the graphene chan- 
nel and nearby defect sites are collected together. If the carriers from 
defect sites can be separated from the total discharge current mea- 
sured at the source/drain, the status of graphene and its environment 
can be systematically investigated. 

Figure lc schematically shows the charge supply and discharge 
process in a graphene channel during a single pulse application. At 
the onset of a negatively biased gate pulse, the potential of the gra- 
phene decreases and hole population in the channel increases (step 
2). This is equal to moving to a higher current point in the current- 
voltage curve of a graphene MOSFET, as shown in Fig. Id. The speed 
of response at step 2 is comparable to the dielectric relaxation time 
around ps. Then, defect sites in the graphene and the nearby inter- 
faces are filled up (step 3). The speed of response at step 2 varies 
depending on the types of defect sites and on the order of lis to ms. 
When the gate pulse is off, the overly populated carriers in the chan- 
nel are discharged through the source and drain contact within a 
dielectric relaxation time (step 4). However, the charges trapped in 
various trap sites are slowly discharged with a time constant in a 
range of us to a few seconds (step 5). As a result, the tail portion of 
the discharge current, 81, highlighted in red in Fig. lc, contains the 
charges released from the defect sites of the graphene as well as other 
surface states near the graphene (to be called charging sites), which 
discharges more slowly than the discharge from the bulk of the 
graphene. 

The baseline current, I 0 , does not change as a function of frequency 
because it is proportional to the total width of the pulse peak as 
shown in Fig. 2a. However, the additional current, 81, discharged 
from the charging sites in each cycle increases as a function of fre- 
quency. Thus, the discharge current can be approximately repre- 
sented using the following equation, 



I c =I„+k-8I-f 



(l) 



where I c is the measured discharge current, I 0 is the current due to the 
charges accumulated in the graphene itself, A: is a frequency- depend- 
ent loss factor to account for the charge loss to the source side due to 
asymmetric metal contacts and other factors, 81 is the additional 
discharge current from charging sites, and / is the measurement 
frequency. 

Typical discharge current measured at 10 KHz to 1 MHz from 
CVD-grown graphene FETs is shown in Fig. 2b. In this work, we used 
large-area monolayer graphene and a number of devices having var- 
iations in the level of physical defect density were fabricated at the 
same time (see Methods and Supplementary Fig. S2). The discharge 
current, I c , is linearly proportional to the frequency in the 10 KHz to 
100 KHz range. As mentioned above, the tail current, 81, consists of 
charges released from defects sites or other surface states. The slope 
slightly increases above 100 KHz and then saturates at a higher 
frequency. 

To understand the frequency response of the I c curve, the charge 
trapping and detrapping mechanism in the graphene should be 
explained first. Recently, our group suggested that the hysteresis of 
a graphene TV curve is primarily caused by two representative 
mechanisms, tunneling and a surface redox reaction 31 . The tunnel- 
ing-induced charge trapping occurs within a few tens of a us time 
constant. The time constant of the surface redox reaction is on the 
order of a few hundred us to ms. Here, the definition of time constant 
is the time to generate 63% of charging. Since the time constant of the 
tunneling process is on the order of 10 us, the tunneling- induced 
discharge current can be modulated up to a few hundred KHz, but 
the discharge current originating from the surface redox reaction 



cannot follow the gate modulation even at a few tens of KHz. 
Thus, the contribution from the surface redox reaction is negligible 
from 10 to 100 KHz. However, this limit does not impede measuring 
defect density because the amount of the surface redox reaction is 
relatively small compared to the tunneling, and it is very weakly 
related to the initial defect density measured using the I(D)/I(G) 
ratio 16 ' 30 . 

In the 10 KHz to 100 KHz range, I c primarily monitors the dis- 
charge current due to the tunneling component, and it is propor- 
tional to the initial defect density as shown below. At this frequency 
range, I c linearly increases in proportion to the frequency. Thus, the 
slope of the linear portion can be used to extract the density of 
charging sites. On the other hand, at a frequency above 100 KHz, 
charge traps that cannot release charges within 5 us start to continue 
holding the charges and those traps no longer contribute to I c . As a 
result, the slope of I c increases and eventually becomes saturated as 
shown in steps 5 and 6 of Fig. 2b. Data obtained from above the 
saturation limit should be ignored because the gate pulse cannot be 
applied correctly to the graphene channel at this frequency due to the 
high channel impedance of the long channel graphene FETs used in 
this work. 

Figure 2c shows I c -f curves for four representative devices with 
different initial defect densities. All four curves show two regions of 
different slopes as a function of frequency at 20 to 100 KHz and 100 
to 150 KHz, respectively. The I c curves are saturated at —200 KHz 
and above. The linear slope at 20 to 100 KHz is attributed to the 
charge exchange between the charge sites and the graphene through 
the tunneling mechanism explained above. To investigate the influ- 
ence of initial defect density (defined as I(D)/I(G) ratio) of graphene, 
we conducted Raman spectroscopy on each graphene channel after 
the graphene channel patterning (Supplementary Fig. S2), and then 
measured the electrical properties after the completion of device 
fabrication. Interestingly, 81 value of thirty five devices are found 
to be closely correlated with the initial defect density, I(D)/I(G) as 
shown in Fig. 2d. This correlation indicates that the 81 value repre- 
senting the density of charging sites originating from both the gra- 
phene bulk channel and various charging sites near the graphene can 
be used to extract the density of graphene defect sites even after 
device fabrication. 

On the other hand, the second slope at 100 KHz to 150 KHz is 
attributed to the change in the base current, I a , due to the Dirac point 
shift during the measurement. Steps 5 and 6 of Figure 2a show that 
the effective charge density at the pulse peak is maintained even 
during the pulse off cycle due to the slow discharge. This is equivalent 
to the parallel shift of TV curves as shown by the blue curve in Fig. Id. 
This gradual increase in the charge density in the channel region 
shifts the Dirac point to the right side and increases I 0 . The second 
slope showed a strong temperature dependence, indicating that it is 
related to the thermally activated mechanism such as chemical 
reaction. (Supplementary Fig. S4) Thus, the second slope should be 
used to monitor the differences in the slow discharge components 
rather than defect monitoring, which represents environmental fac- 
tors such as water molecules trapped around the graphene. 

The slope, 8l/8f, can be translated into density of charging sites 
using following equation, 



^charging site{# / Cm 2 ] = 



2 8l/8f 
k A-q 



(C/t)-t 
cm 2 -C 



(2) 



where A is the device area and N chargjng site is the density of charging 
sites that incurs charging trapping. The value is multiplied by a factor 
of 2 because nearly half the discharge current goes to the source side 
due to a very small drain bias and a symmetrical band structure. Also, 
some of the charges can be dissipated by recombination at the char- 
ging sites. Even though it is not easy to assess the exact amount of 
charge loss, the ratio of charge loss can be simply represented using a 
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Figure 2 | Frequency dependence of discharge current from graphene FETs. (a), Schematic pulse trains showing the amount of discharge current as a 
function of measurement frequency. Step 1 shows total discharge current(I c ) including base current (I„) and discharge current (81) at frequency f. 
As shown in step 1 to step 4, 81 increases as the frequency increases. Step 4 shows the boundary case where the pulse off time is equal to the discharge time. 
Then, the effective pulse on cycle increases effectively and the slow charging mechanism is activated as shown in step 5. Finally, I c is saturated by the 
impedance limit and slow discharge mechanism (step 6). (b), Representative frequency dependence of the discharge current (I c -f) of a graphene FET. 

(c) , I c -f measurements for four representative graphene FETs with different initial defect densities. Inset shows the full range of frequency dependence. 

(d) , The slope of the I c at a frequency lower than 100 kHz (right axis) and the density of charging sites {N c f mrging site ) (left axis) are correlated with the initial 
defect density, I(D)/I(G). 
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Figure 3 | Application of various graphene FETs structures, (a), Device structure (left) and I c -f curves (right) of bottom-gated graphene FETs with and 
without a 30 nm A1 2 0 3 passivation layer. The passivation layer blocks the interaction of graphene with oxygen and water-related adsorbates. (b), Device 
structure (left) and I c -f curves (right) of top-gated graphene FETs on poly(ethylene naphthalate) (PEN) substrates and 90 nm-thick Si0 2 /silicon 
substrates. The 30 nm-thick A1 2 0 3 layer was used as a top dielectric. The graphene on PEN showed much lower defect density than the graphene on 
SiO,/Si. 



constant charge loss factor, k. A typical k factor was 0.6 - 0.7 in 20 to 
100 KHz region for the devices used in this work. 

Finally, Sl/Sf Values are converted to the density of charging sites, 
i.e., ^charging sit^ an d correlated with the initial defect density, 1(D)/ 
1(G) (Fig. 2d), N chargjng site values extracted in this work are on the 
order of 10 14 /cm 2 . This value is relatively high compared to the typ- 
ical interface state density of silicon MOSFETs, which are on the 
order of 10 10 /cm 2 to 10 13 /cm 2 33,34 . However, considering the initial 
quality of CVD graphene is not as good as exfoliated graphene and 
the charging sites include all the defect sites in the graphene channel 
and other trap states near the graphene, this value does not seem 
unreasonable 35 . 

There is appreciable data scattering in Fig. 2d. A part of the data 
scattering is due to the limited accuracy of Raman analysis, which 
represents the quality of graphene only within the beam diameter of 
~ 1 |j.m while the charging current is collected from the whole effec- 
tive active channel area. Also, the scattering in the charging current 
measurements, which is due to the influence of hysteretic device 
characteristics, should be accounted for. Despite these uncertainties, 
the density of the charging sites correlates remarkably well with the 



initial defect density, indicating that the I c -f measurement method 
can be instrumental in analyzing graphene quality, especially after 
device fabrication, because this method can be used to compare the 
impact of processes on the same graphene even without calibration. 

The DCA method may look similar to the charge pumping 
method, which measures the interface state density of silicon 
MOSFETs, but the measurement principle is quite different 31,32 . In 
the DCA method, the slope of the discharge current- frequency (I c -f) 
curve is used to analyze the density of the charging sites while the 
electron-hole pair recombination at the interfacial charging sites is 
measured at a fixed frequency in the charge pumping method. 

Figure 3 shows examples of various graphene FETs structures 
analyzed using the DCA method. Figure 3a shows the I c -f curves 
measured in bottom gate graphene FETs with/without a 30 nm 
A1 2 0 3 passivation layer. Since the A1 2 0 3 passivation decreases con- 
tamination from the graphene surface, the density of the charging 
sites decreased to 2.47 X 10 14 /cm 3 from 3.19 X 10 14 /cm 3 after pas- 
sivation. A more prominent slope increase before the passivation in 
the 100 to 180 KHz region indicates that the graphene surface is 
more prone to the water redox-induced charge generation, as 
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expected 25,36 38 . A second example compares the top gate graphene 
FETs on a 90 nm Si0 2 /silicon substrate with a poly( ethylene 
naphthalate) (PEN) substrate (Fig. 3b). In this case, 30 nm A1 2 0 3 
is used as the gate dielectric layer for both device structures. Since the 
PEN substrate is hydrophobic, the I 0 shift due to the water redox 
reaction is expected to be much less pronounced than that of the Si0 2 
substrate 8 . Also, more stable I-V characteristics have been reported 
for hydrophobic substrates 39 " 41 . Thus, the density of charging sites 
should be lower with graphene on PEN substrates. Even though the 
initial quality of the graphene was similar, the density of charging 
sites in the graphene FET on PEN —6.25 X 10 13 /cm 2 was much less 
than —6.66 X 10 14 /cm 2 of the Si0 2 /Si substrate. Even though the 
charge loss factor, k> was very low —0.2 for top gate graphene FETs 
because of fast recombination in the charging sites at the A1 2 0 3 
interface, a relative comparison between different substrates was 
feasible. The much lower N charsit1g site in the PEN substrate is prim- 
arily attributed to the reduced interaction between the graphene and 
PEN substrate. These examples clearly demonstrate that the new 
DCA method can be used to assess the influence of environmental 
factors and the quality of graphene after device fabrication. 

Even though this unique electrical characterization method is the 
first attempt to quantitatively analyze the density of charging sites of 
graphene FETs after device fabrication, the preceding examples show 
how instrumental this method can be in optimizing graphene device 
structures. Similar to the concept of charge pumping for silicon 
MOSFETs, which has evolved into several advanced CP methods, 
there might be several modifications that can be made to improve the 
accuracy of this method and enhance the quality of information 
extracted from graphene FETs. Yet, the DCA method proposed in 
this work will provide an important stepping stone towards the real- 
ization of high quality graphene FETs. 

Methods 

Methods summary. Experimental graphene samples were prepared through transfer 
process and conventional photolithographic technique using graphene sheets grown 
on copper foil by a CVD method. Before the completion of the device fabrication, 
Raman spectroscopy was performed on the patterned graphene channel. Electrical 
properties of fabricated graphene FETs are examined with continuous pulse system. 
In order to monitor the defect density of graphene in a complete device structure, 
discharge current was measured by applying square pulse biases at 10-1,000 KHz 
using a pulse generator and a parameter analyzer. 

Graphene synthesis and transfer. A 1 cm X 1 cm monolayer of graphene sheet 
grown on Cu foil using a chemical vapor deposition (CVD) process was transferred to 
90 nm Si0 2 thermally grown on P-doped Si substrate using a PMMA- mediated 
transfer method 29,30 . The optical photograph and the Raman data of graphene used in 
this work showed that the graphene is mostly a monolayer with the low defect density 
(Supplementary Fig. S2). 

Fabrication of graphene FET devices. For the device fabrication, 100 nm Au 
electrode patterns were formed on a graphene sheet using i-line contact 
photolithography and a lift off process. Then, graphene channels between source and 
drain electrode were patterned using i-line photolithography and oxygen plasma 
ashing process (process power = 50 W, process time = 90 seconds). The surface of 
graphene channel was passivated with 30 nm A1 2 0 3 using atomic layer deposition 
(ALD) process at 130 C and annealed in N 2 ambient at 200 C for 30 min to minimize 
the influence of ambient. 

Raman spectroscopy and mapping. After the graphene channel patterning, the 
quality of graphene channel defined as I(D)/I(G) ratio was measured for all devices 
using Raman spectroscopy (Renishaw, X — 248 nm, power = 20 mW). 
Supplementary Fig. S2b shows the representative Raman spectra of the single-layer 
CVD graphene. To investigate the influence of the initial defect density, 35 graphene 
channels with various I(D)/I(G) ratio from 0.13 to 0.41 are chosen from various 
locations of graphene sheet as shown in Supplementary Fig. S2c and S2d. Electrical 
characteristics of these devices are monitored after the completion of device 
fabrication. (Supplementary Fig. S3). 

Discharge current measurements. For discharge current analysis, square pulse 
biases in 1 0 KHz to 500 KHz range were applied to the substrate, which works as a 
bottom gate electrode with 90 nm Si0 2 gate dielectric, using a pulse generator 
(Agilent 81110) as shown in Fig. la. Pulse width was determined by the duty cycle of 



pulse ( — 50%) at a given frequency and the rise and fall time was fixed as 100 ns. The 
source is connected to a ground and the drain is connected to a parameter analyzer 
(Keithley-4200 SCS) for a current measurement. 
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